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Abstract. We investigate the scattering of 2D cylindrical invisibility cloaks with 
simplified constitutive parameters with the assistance of scattering coefficients. We 
show that the scattering of the cloaks originates not only from the boundary conditions 
but also from the spatial variation of the component of permittivity /permeability. 
According to our formulation, we propose some restrictions to the invisibility cloak 
in order to minimize its scattering after the simplification has taken place. With 
our theoretical analysis, it is possible to design a simplified cloak by using some 
peculiar composites like photonic crystals (PCs) which mimic an effective refractive 
index landscape rather than offering effective constitutives, meanwhile canceling the 
scattering from the inner and outer boundaries. 
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1. Introduction 

Keeping the form of Maxwell equations invariant, coordinate transformation (CT) 
provides us the possibility to control the electromagnetic field, by changing the 
distribution of permittivity and permeability, as well as to create an isolated space 
to hide any objects, i.e. invisibility cloaking [Tl|2]. In recent years, invisibility cloaking 
has been attracting significant attention [31 HJ [5l El [71 [H HI [TOl [TTl [12] , due to its novel 
electromagnetic (EM) characteristics and potential applications. 

In general, the invisibility cloaks, realized through a CT approach, possess 
constitutive parameters being anisotropic and spatially varying. In particular, a 3D 
spherical invisibility cloak calls for constitutive parameters which are anisotropic and 
inhomogeneous, but with finite values, which seems realizable in practice [1]. However, 
in the 2D case, cylindrical invisibility cloaks commonly possess constitutives not only 
being anisotropic and spatially dispersive, but also requesting infinite components [1], 
which in practice leaves the realization of any ideal cloak impossible. To overcome this 
drawback, Schurig et al. suggested the simplification of the constitutive parameters, i.e. 
flexibly choosing the constitutives while keeping the spatial distribution of dispersion 
relations unchanged [1], which is based on the ray tracing concept. With such an 
approach, the first cylindrical cloak shell was successfully constructed. 

Since we have to simplify the constitutives of a 2D cylindrical cloak in the practical 
design and fabrication, the scattering of the cloak itself is unavoidable, i.e. an invisibility 
cloak with reduced parameters is inherently visible [13|. In order to reduce the scattering 
of a cloak with reduced parameters, a variety of research work has been done. For 
instance, a second-order polynomial function was used in the transformation in order 
to avoid the reflection at the outer boundary [H], or to minimize the scattering of 
cloaks with reduced parameters by the optimized simplification [151 (HI (Uj- Though 
people have achieved success in the invisibility cloak theory and experiments, a general 
description of the scattering by a cloak with reduced parameters is unrevealed. 

In this paper, we analytically solve the scattering problem of a cylindrical invisibility 
cloak by applying the classical scattering theorem. With the aid of scattering coefficients, 
we show the full-wave solution of the scattering by cloaks with ideal as well as reduced 
parameters. According to our theoretical analysis, the scattering of a cloak with reduced 
parameters is contributed by both the inner and outer boundary, as well as the spatially 
varying internal impedance. We point out the way to minimize the scattering of a 
cloak in practice, in which case its scattering is dominated by the spatially distributed 
permittivity/permeability rather than the boundary conditions. In addition, the way to 
realize nearly zero scattering by an approximate cloaking shell is also proposed. 

2. Scattering of cylindrical invisibility cloak 

Consider the anisotropic and inhomogeneous invisibility cloaks obtained from a CT 
approach. If the mapping employed in the transformation is simple, i.e. only one 
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coordinate variable is involved in the transformation like the ones proposed in Ref. [T] 
or Ref. |1] , the cloak's e and /2 are diagonal and inhomogeneous. Though the analytical 
solution of the scattering problem of an ideal cloak can be found by applying the Maxwell 
equations, the scattered fields by those cloaks with reduced parameters have to be 
determined numerically pSj. However, the scattering by a multi-layered cloak shell can 
be found by applying the full wave expansion, satisfying the boundary conditions |18j . 

Suppose there is an A^-layered cloak structure as shown in Fig. [1] (A). The inner 
and outer radius of this cloak shell are a and 6, respectively, pj is the outer radius of 
the j-th shell, po = a stands for the radius of the inner region, pj^ = b is the outer 
radius. Inside the j-th layer with thickness dj = pj — see Fig. [H the permittivity 
and the permeability are constant and denoted by ej and Pj . If ej and pj vary 
continuously in space, then a continuous cloak shell can be obtained if we let iV — )• oo 
and maxdj — )■ (j = l,2,3...iV ). For simplicity, the permittivity and permeability in 
each layer of the multi-layered cloak are assumed to be anisotropic but diagonal, e.g. 
ej/eo = diag e^j, e^j, j . Here, we consider the TM polarization case, i.e. the electric 
field is polarized along z direction, while the same procedure can be applied to a TE 
polarization case. 

Basically, the field inside each layer can be easily expressed by a superposition of 
vector cylindrical waves [191 H] • We rewrite the field inside the j-th layer as 

oo 

E.,j= E [CirnJu,iK,p)+Ci^Hj'\k,,p)]e^"^'^, (1) 



Figure 1. Configuration (A) and Cross section (B) of a mufti-layered cyfindricaf cfoafc. 
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h ■ 



and fco = w^eoA^o- Here, 



^ and 6*2 ^ are the unknown expansion coefficients inside the j-th layer which need 



to be determined by matching the boundary conditions. In the outer space, Cf^^^ and 
C2,m^ are used to represent the total field in the outer free space, composed of the 
incident and the scattered fields. In Eq. ([1]), the first term represents the standing wave, 
while the second one is the scattered wave in each layer |T9j. Boundary conditions are 
applied to determine the unknown coefficients. At the interface between the j-th and 
the j + 1-th layers, tangential electric and magnetic fields should be continuous, hence 
we find 

^ l,m 



with 



mj+1 



m,j+l^ m,j 



Qr. 



a 



2,m 



^Hii;{k,,^,p,] 



<t>,3 



Q 



•JujikpjPj) 

'J'uX^pjPj) 



'--Hl'}\k 



(2) 



(3) 



(4) 



Similarly, we consider the boundary conditions at any other interface, and the two 
unknown coefficients in the free space can be evaluated by 



Pm,N+l 
Qm,N+l 



^ --1 

UP 

7 = 1 



+l^m,j 



Qr 



^1 



(5) 



Equation is the common relationship that a multi-layered cloak should satisfy, 
with which the scattered field of a cloak consisting of several layers can be conveniently 
worked out [IB]. However, since the 2D invisibility cloak requires that k^f^j tends to 
zero near the inner interface, then uj is singular which makes both Ti ^ and T2,m being 
ill-conditioned and the numerical evaluation of Eq. (|5]) fails. 

However, since our purpose is to find out the scattered field, the scattering 
coefficients are of course very helpful and simplifying the formulation may also overcome 
the drawback of Eq. (jS]) . The scattering coefficient for the m-th cylindrical wave in the 
j-th layer is defined as Rl^ = Cl^^ICi^^^. From Eq. (|2]), we easily arrive at 

where A^, 5^, C^, and Z)^ are functions of z/, fcp, k^ in both the j-th and j + 1-th 
layers. If we further assume that all the layers are very thin, and let — )■ oo and 
maxfij — )■ 0, i. e. the cloak structure gradually tends to be continuous, we can rewrite 
Rl^\ Ri, Al, Bl, CI and as + A/?„, R^, A^ + AA^, + AB^, C^ + ^C^ 



and Dm + ^Dm, respectively. Hence we find 



ARr 



Am + AAm + {Bm Cm + ABm ACm)Rm {Dm + ADm)Rn 
Cm + + {Dm + ADm)Rm 



(7) 
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Balancing the first order differentials in Eq. ([7]), we get 



(8) 



with 



kpk'ppJ'l{kpp) + k'pJlikpp) + kpv'—Jl{kpp) 



+J'u{kpp) 



kpk'ppJlikpp) + kp^J^{kpp) + kpiy'-^J^{kpp) 



B'^ = -Hl'\kpp) 



d 



kpk'ppj'l{kpp) + k'pJl{kpp) + kpu'—Jl{kpp) 



+Hi'Hkpp) 



kpk'ppJlikpp) + kp—Jy{kpp) + kpp'—Jy{kpp) 



C'm — Jv{kpp) 



kpk'ppHi'y\kpp) + k'pHi'y{kpp) + kpu'—H^^y{kpp) 

-Jlikpp) kpk'ppH^y {kpp) + kp^Hi'\kpp) + kpv'^H^\kpp) 



D'^ = Hl'\kpp) 



kpk'ppHi'y\kpp) + k'pHi'y{kpp) + kpv'^H^;y{kpp) 



-Hi^yikpp) Ikpk'ppHl'yikpp) + kp^Hl'\kpp) + kpu'^Hl'\kpp) 

For more details of the technical work, please refer to the appendix. Containing the 
term involving Eq. (j8]) becomes the non-linear Riccati equation for continuous 

invisibility cloaks. It could be solved numerically with initial value at the boundaries. 

To find out the scattering coefficients in the outer free space, we should find 
out the scattering coefficients on the inner interface first. For the ideal cloaks, the 
permittivity/permeability on the inner boundary is infinite. This prevents us from the 
evaluation of the scattering coefficients, hence we have to assume that the cloak's inner 
radius has a small perturbation 5, i. e. the inner radius is now a + 5 but the rest of the 
parameters are kept unchanged. Notice that there is no scattered field inside the inner 
region (p < a + 5), then the scattering coefficients on the new inner boundary can be 
evaluated by applying Eq. (|2]) 

n , , Va+5Jm[ki{a + 5)]Jl^^^[ka+5,p{a + 5)] - Ji[ki{a + S)]J^^^^[ka+5,p{a + 6)] 

Hm{a+d) = -— : r: — r— n^ : tt' y^) 



Va+5Jm[ki{a + 6)]H!)'J_^, [ka+s,p{(^ + ^)] - JL[h{a + 6)]H!)^X,[ka+5,p{(^ + ^)] 



with r)a+5 = \J ^a+5,zl Pi'a+s,(t> and ki being the wave number in the region p < a + 6. 
The initial scattering coefficients of any cloak structure can be obtained by letting 
(5 — )■ in Eq. ([9]). For example, for ideal cloaks, — )■ and Va+s — )■ lead to 

Ma + 6)^ -T#il7r^ and R^a + 6)^ 0(m ^ 0). 

To verify our theoretical formulation, here we assume there is a cylindrical 
invisibility cloak with a=0.024m, 6=0. 072m and its constitutive parameters = pi^ = 
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diag (^)^(^-^) • By using Eqs. ([HD and ([21), the scattering coefficients can be 

conveniently evaluated. Fig. [2] (A) shows the scattering coefficients for the lowest four 
cylindrical waves (m=0,l,2,3). 

From Fig. [21 we see that all the scattering coefficients (|-Rm|) are zero on the outer 
boundary. Since an ideal cloak perfectly matches the outer space, no scattered wave 
will be excited. Likewise, the inner boundary, \Rm\ is also zero, except the 0-th order 
case. Though non-zero no electric field can be excited inside the inner region [1], 
which can also be proved by applying Eq. ([21). In Fig. [21 (A), Rq varies fast near the 
inner boundary, hence the total scattered field is sensitive to any perturbation of the 
inner boundary, as reported in Ref. [20] . 



3. Cloaks with simplified parameters 

In the previous section, we have obtained the scattering coefficients for an arbitrary 
invisibility cloak structure with cylindrical symmetry. However, the permittivity and 
permeability of ideal cloaks are not only anisotropic and inhomogeneous, but also have 
infinite components (e^ and /i^) on the inner boundary j3]. Hence, in order to realize 
the practical design and fabrication, we have to simplify the constitutive parameters, i. 
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Figure 2. (Color online) (A) The absolute scattering coefHcients for an ideal cloak 
with linear transformation function. (B) The total near field distribution with a bare 
perfect electric conducting (PEC) cylinder scattering. (C) The near field distribution 
with a PEC cylinder covered by an ideal cloak. In the calculation, operation frequency 
is 7GHz. 
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e. maintain the dispersion relation distribution, while releasing our choice of e and /i 
to a more realistic set [H |2T1 |T5] . Unfortunately, this simplification generally affects 
the scattering coefficients, i.e. the mismatched inner and outer boundaries, as well as 
the variation of (for TM case), will cause non-zero scattering coefficients, which we 
see below in detail. 

In Ref. [H], to minimize the scattered field after the simplification, high order 
coordinate transformation was applied, for the purpose of impedance matching on 
the boundary p = b. While the inner boundary is usually assumed to be perfectly 
conducting O [HI [18] , by doing this the inner region is isolated from the outer space. 
However, the 0-th order cylindrical wave can still see the inner region and detect the 
perfect conducting wall [13]. Although with the help of some numerical method like 
genetic algorithm (GA) we can optimize the scattered field by a multi-layered cloak [TH] , 
the physical interpretation of the scattering is missing. Returning to Eq. ([9]), the initial 
value of R is commonly determined by both the cloak and the object placed inside. 
However, we notice that R{a) can be independent of the permeability and keeps zero 
if fca,p = 0, which means that the scattering of a cloak will be irrelevant to the inner 
objects. Again, in Eq. (|8]), the alternative /x^ after the simplification may introduce the 
disturbance of scattering coefficients through the term containing /i'^/ /i^. To minimize 
the contribution from this term, the other terms should dominate, i.e. k'^ can not be 
zero. 

If the coordinate transformation is p' = f{p), (p' = (p and z' = z, all the above 
restrictions can be expressed as f'{b) = 1, /'(a) = 0, f"{p) ^ 0, with natural 
conditions /(a) = and f{b) = b. Any functions satisfying these restrictions can 
be applied to the cloak design. Here, a third-order polynomial function is adopted, i.e. 
p' = Ap^+Bp^+Cp+D with A = -(^, B = 2(feb) -^^A C = l-3Ab^-2bB, and 
D = —{Aa^+Ba'^ + Ca). We should notice that the choice for higher-ordered polynomial 
function is not unique, which indicates that the optimization might be accomplished by 
using some numerical method like GA. 

Now we address some numerical calculations to verify our theoretical analysis. We 
suppose an infinitely long cylinder with radius r=24mm made of PEC is placed in 
free space. The scattering coefficients of this PEC cylinder under TM polarized wave 
illuminating are listed in the second column of Tab. [H while the near field distribution 
is shown in Fig. [2] (B). 

Next we assume that a cylindrical cloak with a=24mm and 6=72mm is applied 
to hide the specified PEC cylinder. First, we apply a linear function to complete the 
CT, as previously stated. In practice, the constitutive parameters have to be simplified. 
Here we make p^=l while we tune the other parameters correspondingly |1]. The near 
field distribution of the scattering by this simplified cloak is shown in Fig. [3] (A), while 
the scattering coefficients for the lowest four order cylindrical waves are listed in the 
third column in Tab. [TJ We clearly find that the non-zero scattering coefficients make 
the cloak visible [I3]. Comparing the second and the third columns in Tab. [H we see 
that though with non-zero scattering coefficients, this simplified cloak has a larger size 
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Table 1. The absolute scattering coefficients for tire lowest four cylindrical waves are 
shown. In the calculation, /=7GHz. \Rc\'- the bare PEC cylinder case; \Ri\: PEC 
cylinder covered by a simplified cloak from linear CT; |i?3or(ih PEC cylinder covered 
by a simplified cloak from a third-order polynomial CT; |i?apph PEC cylinder covered 
by an approximate cloak shell. 



m 


\Rc\ 


\Ri\ 


\R3ord\ 


1 Rap'p 1 





0.9036 


0.801 


0.354 


0.1187 


1 


0.3004 


0.197 


0.135 


0.0068 


2 


0.9934 


0.082 


0.031 


0.0066 


3 


0.7418 


0.245 


0.060 


0.0120 



but less scattering than the PEC cyhnder. 

Second, a third-order polynomial function, which can be uniquely determined by 
the above analysis, is applied in the CT. Though the cloak still needs simplification, its 
scattering is only caused by the alternated /x^, see Eq. (IHl), hence we would expect the 
reduction of all the scattering coefficients. We still make after the simplification. 
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Figure 3. (Color online) Total electric field distribution for a PEC cylinder covered 
by (A) a simplified cloak coming from a linear CT; (B) a simplified cloak coming fi'om 
a third-order polynomial CT; (C) an approximate cloak shell. (D) The magnitude of 
the scattered electric field on the virtue contour C (p=0.084m). Due to the symmetry, 
here we only show the distribution in < (/) < tt. In all of the calculation, /=7GIIz. 
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The total field distribution is shown in Fig. |3] (B), while the lowest four scattering 
coefficients are listed in the fourth column in Tab. [H We see that though the scattering 
coefficients are non-zero again, this cloak scatters less EM power than the one shown in 
Fig. [3] (A). 

From Fig. [3] and Tab. [H a simplified cloak structure related to a third-order 
polynomial CT has less scattering than the one related to a linear CT. However 
the scattering contributed from the /i^ term in Eq. ([8]) is still unavoidable. In 
fact, by applying high-order polynomial transformation functions and some numerical 
optimization, the cloak's parameters can be optimized to minimize the scattering 
involved by the /i<^ term. Here, we point out one possible analytical solution. 

We go to Eq. ([8]) and if the term containing fi'^/ fi^f, is zero, then the scattering of 
the cloak is dominated by the spatially dispersive refractive index. For this purpose, 
we have ^'^/ fiip = 0, or p' = b^~'^p^ , with X being determined by boundary conditions. 
Rigorously speaking, functions of this kind can not become the transformation function 
since p' 7^ if p 7^ 0. However, properly choosing the two unknown coefficients can 
still support an approximate cloak shell [22j. An approximate cloak shell coming 
from this kind of transformation function will possess constitutives as = jf = 
diag X~^, X, Xfe^^^^p^*^"^ , i.e. no singularity existing in both the permittivity and 
the permeability. Since this kind of invisibility cloak shell has been reported in Ref . [22] , 
we do not place too much discussion here. The scattering coefficients are summarized in 
the last column in Tab. [1] and the electric field distribution is given in Fig. [3] (C). We see 
that the approximate cloak shell scatters little EM power, which in principle "shrinks" 
the non-zero scatter [221. 



4. Conclusion 



In conclusion, we have solved the scattering problem of 2D cloak structures with 
cylindrical symmetry, for both the ideal and simplified cases. We show that the 
scattering coefficients of a cloak are determined by both the inhomogeneous refractive 
index and the inhomogeneous p^(e</,), as well as the inner and outer boundary conditions. 
In the practical fabrication, to avoid the infinities of an ideal cloak, the cloak's 
constitutive parameters have to be simplified, which commonly makes an otherwise 
invisible cloak visible. 

We point out the way to avoid the scattering from both the outer and the inner 
boundaries, i.e. to have the outer boundary matched to the free space and make 
the refractive index tend to zero to minimize the scattering from the inner boundary, 
e.g. high-order polynomial CT can be applied to the cloak design. However, though 
numerical results show the reduction of scattering, the third-order polynomial based CT 
can not guarantee the zero scattering due to the changing of Ptf> after the simplification. 
Some numerical method like GA can be adopted to further optimize the scattering of 
a simplified cloak. Our theoretical work shown in this paper clearly exhibits the role 
of the inhomogeneous refractive index. With minimizing the scattering involved by the 
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jj,^ term and the two boundaries, dielectric metamaterials and photonic crystals, which 
prefer effective spatial dispersion rather than effective constitutive parameters, could 
be conveniently applied to construct an invisibility cloak [231 Ell ES], in which case 
the topology optimization may stimulate the practical design and fabrication |26j. Our 
theoretical work in this paper would improve the practical fabrication and application 
of invisibility cloaks in future. 
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Appendix A. The derivation of Eq. 



We rewrite Eq. ^ as 
and after some algebra operation we find 



TT 

2i 



m m 
in m 



(Al) 



^in = ~^p,i+iPjJvj+ii.^p,j+iPi)Juj{^PdPi) + ^p,jPj~~? J^j+i{^p,i+^Pi)Juj{^p,jPj) 



P<p. 



Bin — ^pj+iPjJvj+ii.^p,j+iPj)H^^ {kpjPj) + kpjPj—^ Jy.^^{kpjj^iPj)Hl}^ (kpjPj), 

P<t>,j 

cL = ^p,i+iPi^^J|/(^p,j+iPj)^(^pjPj) - K,jpj^7^Hl]l^{kp^^+ip^)Ji^{kpjpj), 

p<f>,j 

DL = K.+iP^H^l'ikp.+iPM]\K,P^) - fcp.P.^<!.(fcp.+iP.)<^'(^P.P.)- 

p<l>,j 

From Eq. (lA.ll) . we can derive Eq. ([6]). Furthermore, assuming that all the cylindrical 
shells are thin, i.e. max dj — )■ 0, we can rewrite R{^^, Rln, Al^, i?^, and Z)^ as 
Rm + ARm, Rm, + AA^, 5„ + A5„, Cm + AC^n and Dm + ADm- The A sign 
appears in those terms involving Ap = pj^i — pj. Thus Eq. (El) has a new form as Eq. ([71) 
which is rewritten here 

Am + AAm + {Bm — Cm + ABm ~ ACm)Rm ~ {Dm + ADm)Rm 



AR„ 
with 



Cm + ACm + {Dm + ADm)Rn 



(A.2) 



A 

B. 



m 0, 

2i 

IT ' 
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Dm = 0. (A.3) 
In deriving Eq. (IA.3P , we use the Wronskian of Hankel functions [27j 

Hi'\xMx) - Mx)Hi'H^) = • (A.4) 

TTX 

Furthermore, simphfying Eq. flA.2p . we get 

ARm{Cm + ACm + ADmRm) = AAm + {ABm- ACm) Rm- ^DmRm^ .{A.5) 

Balancing the first order differentials in Eq. (lA.Sp . we can get Eq. ([8]). 
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